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ABSTRACT 

A search for matched circle pairs of similar temperature fluctuations in the final 
WMAP 9yr data is carried out. Such a signature is expected if the space of the 
Universe is multi-connected. We investigate the relation between the pixel resolution 
of CMB maps and a Gaussian smoothing in order to lower the probability for missing 
matched circle pairs. CMB maps having the 3-torus topology are generated with the 
characteristics of the WMAP satellite in order to determine how large the smoothing 
should be chosen in CMB maps disturbed by detector noise. The V and W band data 
are analysed with respect to matched circle pairs and a tentatively signal is found 
for a circle pair, which lies, however, close to the plane of the Galaxy. It is, however, 
inconclusive whether this signal is generated by chance, is due to residual foregrounds 
contained in the V and W band maps, or is due to a genuine topology. 
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1 INTRODUCTION 

Since the cosmic microwave background (CMB) radiation 
provides us with the earliest admissible electromagnetic 
radiation, it contains information about the largest scales 
of our Universe. Therefore, it is most promising to search 
for topological signatures in the CMB radiation in order to 
reveal the topology of our Universe. One of t he topological 
tests is the matched circle test proposed by I Cornish et al.l 
(1998) on which this paper puts its focus. For an intro- 
duction into the topic of cosmi c topology and discussions 
concer ni ng topological tests, see iLachieze-Rev and Luminetl 



19951); lEuminet and Roukemal (Il999l); iLevinl ||2002|) 



iReboucas and Gomerd J2004l ; Luminetl l|200Sl ); Mota et al.l 



( 201(1 l201ll ); iFuiii and Yoshiil pOllD . 



The idea behind this test is the following. The CMB 
sky observed from a given observer position originates from a 
sphere around this observer, i. e. from the surface of last scat- 
tering (SLS). If the Universe possesses a non-trivial topology, 
the space can be viewed as being tessellated by cells which 
have to be identified. This space is the so-called universal 
cover. Each cell has a "clone" of the observer, which ob- 
serves the same CMB sky. If the observer and a given clone 
are not farther separated than the diameter of the surface 
of last scattering, the two spheres overlap in the universal 
cover. The intersection of the two spheres is a circle seen by 
the observer and its clone in a different direction. Since the 
observer and its clone are to be identified, one concludes that 
two circles should exist on the CMB sky with identical tem- 
perature fluctuations seen in different directions but with 



the same radius. A non-trivial topology is thus betrayed by 
as much pairs of circles with the same temperature fluctu- 
ations as there are clones of the observer not farther away 
than the diameter of the surface of last scattering. 

This description is simplified since it ignores several ef- 
fects which alter the CMB temperature fluctuations of the 
observer and its clone in a different way, so that the tem- 
perature fluctuations on both circles are no longer identi- 
cal. The two most important CMB contributions in this 
respect are the Doppler contribution, whose magnitude 
depends on the velocity projection towards the observer, 
and the integrated Sachs- Wolfe (ISW) contribution, which 
arises along the path from the surface of last scattering 
to the observer or to the clone. These two paths are not 
identified and lead to different contributions to the to- 
tal CMB signal. These i ssues are discus s ed in more detail 
by I Aurich et al.l <|2006|); iRiazuelo et al.l i|2006l ); iKev et all 
|2007l) ; iBielewicz and Bandavl (|201ll v" There are further de- 
grading effects, e. g. the finite thickness of the surface of last 
scattering, but these are considered as subdominant to the 
Doppler and ISW contribution. Despite these degrading ef- 
fects, there should be a detectable topological signal. 

In addition to these contributions, which can be com- 
puted with the standard CMB physics, there are residuals 
left over by the subtraction of foreground sources which have 
their own uncertainties. 

The circles-in-the-sky (CITS) test requires a full sky 
survey and has been applied t o different sky map s de- 
rived from the W MAP mission. I Cornish et all (|2004l ) and 
IKev et all (|2007r i analyse the first year CMB data with re- 
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spect to nearly back-to-back circle pairs with a negative re- 
sul t. A hint of the P oinca re dodecahedral topology is found 
bv iRoukema eTall (Efflol l in the first year data. The cir- 
cle signals correspond to circles with a sm all radius o f abou t 
11°. This claim could not be confirm ed bv lKev et all (|2007l ). 
The analysis of I Aurich et~ai1 (|2006l ) uses a weight function 
constructed to suppress the Doppler and ISW contributions. 
It does not found a convincing signal, and only a marginal 
hint is found in favour of the Poincare dodecahedral space. 
The three year WMAP data are inves t igated with respect to 
the 3-torus topology bv I Aurich et al.l l|2008j). and it is found 
that the degree of uncertainties in the CMB map must be 
significantly below 50/iK in order to have a realistic chance 
to discover such a topology. The seven year WMAP data 
are analysed bv lBielewicz and Bandavl l|201ll ) with respect 
to back-to-back circle pairs. They rule out topologies having 
such circle pairs with radii larger than 10°. This work em- 
phasises the presence of residual Galactic foreground emis- 
sion close to the Galactic plane in the WMAP ILC 7yr map 
and chooses to use a mask to eliminate this non-CMB con- 
tribution. The search for gener al circle pairs that are no t 
back-to-back is carried out by IVaudrevange et"al] l|2012l ). 
Although some signals are detected, they are ascribed to 
the residual foreground contained in the ILC map which is 
used within the WMAP Kpl2 sky mask. This analysis con- 
cludes that no topological signal is present in the seven year 
WMAP data with circle radii above 10° . A CITS search for 
an or bifold line topology is carried out by iRathaus et al.l 
i|2013h and a tentative candidate match is found. 

Are those topological spaces definitely ruled out by 
these CMB analyses? In order to exclude or find possible 
loop holes in the matched circle pair signal, we generate a 
high resolution CMB sky map for the 3-torus topology which 
requires a space with zero curvature, i. e. an Euclidean space. 
This map is analysed with respect to the amplitude of the 
topological signal. 

The 3-torus map is computed using the cosmological 
parameters of the ACDM conco rdance model which are 
published bv lJarosik et al.l (|201ll ) in their Table 8, column 
"WMAP+BAO+i/o" (ttbar = 0.0456, Q. cdm = 0.227, Hub- 
ble constant h = 0.704, reionization optical depth r = 0.087, 
scalar spectral index n s = 0.963). In the following only 
the cubic 3-torus is considered where all three topologi- 
cal lengths have the same value L. In units of the Hubble 
length Lh = c/Ho, the chosen side length is L = 1.5 in 
sections [2] to [4] This value is significantly smaller than that 
obtained from the requirement that the CMB temperature 
correlations of the model should match the observed corre- 
lations. A best fit of th e correlations leads to val ues around 
L — 4, see figure 11 in lAurich and Lustid (|201ll ). It should 
be noted that this 3-torus topology with L — 4 could be de- 
tec ted in future high-reds hift galaxy surveys as emphasised 
bv IRoukema et al.l (|2013h which would be independent of 
the CMB analyses. The reason for the unrealistically small 
value of L = 1.5 is that the CMB anisotropics have to be 
computed by using the eigenmodes of the 3-torus, and we 
want to accurately compute the CMB up to the multipolc 
moment Ci with i max = 3000. To achieve this fine-structure 
for the chosen cosmological parameters, all eigenmodes with 
a wavenumber k up to k — 1 026 are included. This means 
that 61 556 892 different wavenumbers k are taken into ac- 
count. The simulation of a CMB map for L = 4 with that 



resolution is currently not feasible. However, the amplitude 
of the signal due to a single pair of matched circles should be 
independent of the size of the fundamental cell. The size de- 
termines, however, the number of matched circle pairs which 
decreases with increasing size of the 3-torus cell. The cubic 
torus topology with L — 4 possesses six circle pairs with 
radius a ~ 31° and three circle pairs with radius a ~ 53°. 
The total number of circle pairs increases for L = 1.5 to 182. 
Thus, the simulation with L = 1.5 has the advantage that 
there are much more circle pairs to analyse with respect to 
their discovery probability. For the decision, whether a cell 
with L = 4 is discovered with certainty, one has to take the 
much smaller number of circle pairs into account. The CMB 
sky map for the 3 -torus with L = 1 .5 is generated in the 
HEALPix format (|G6rski et all l2005f ) with a resolution of 
A^idc = 4096. This map possesses a pixel size of 51". 5 which 
is three times smaller than the resolution of the simulation 
being of the order 180°// ma x ~ 3'. 6. The following analysis 
is mainly based on this map. We will also use less accurate 
L = 4 simulations in sections [S] and |S] when the detector 
noise is taken into account since in that case a lower accu- 
racy of the simulation is sufficient. Furthermore, we restrict 
us to a search of back-to-back circle pairs. 



2 THE CITS AMPLITUDE FOR INDIVIDUAL 
CIRCLE PAIRS 

Let us now turn to the correlation measure used for the 
detection of CITS signals. The expansion of the temper- 
ature STi(<j)) along a circle in a Fourier series 5Ti(<f)) — 



< 2-7T, allows one to define the m- 



weighted cir cle signature for tw o circles i and j having a 
radius a as i|Cornish et al.ll2004h 



2E 



,im/3 



y 
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The angle /3 takes a possible shift between the two circles 
into account. A perfect correlation would be revealed by 
Sij(a,/3) = 1 but due to the degrading Doppler and ISW 
contributions, the actual value is lower. Then the maximum 
is taken over all circle pairs and all shift angles j3 for a fixed 
circle radius a, i. e. 



5(a) 



max Sij(a, (3) 

i,j,0 



(2) 



In the following we restrict this analysis to back-to-back 
circles which determines the index j as a function of i. Every 
pixel i of the HEALPix map is treated as a possible centre 
of a circle in the CITS search. 

Even if the correlation S(a) is computed for a CMB 
map belonging to a model with a simply-connected topol- 
ogy, one obtains non-vanishing correlations by chance. The 
smaller the rings, the higher is the probability that the tem- 
perature fluctuations are accidentally correlated. This leads 
to a background for S(a), which starts at S(a) ~ 1 for a — 
and declines towards larger radii a, until at a — 90° the pairs 
are mapped onto each other leading again to S(a) ~ 1. A 
topology can only be discerned if the CITS correlation due 
to a matched circle pair is significantly larger than this back- 
ground. 

This background behaviour is shown in figure [TJ where 
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Figure 1. The CITS correlation 5(a) is calculated for a 3-torus 
simulation with L = 1.5 using the HEALPix resolution N s ^ e = 
512 and a Gaussian smoothing with 8 = 0.5°. The dotted curve 
shows S(a) obtained from all back-to-back circle pairs and all 
shift angles j3. The open circles reveal the CITS values obtained 
from the 182 matched circle pairs of this topology. The dashed, 
full, and dotted-dashed curves belong to the subset of circle pairs 
having the same circle centres as one of the matched circle pairs 
with radii a = 20.5°, 24.7°, and 56.2°, respectively. The CITS 
values of these matched circle pairs are depicted by red dots. 



the CITS correlation S(a) is computed from a 3-torus sim- 
ulation with L = 1.5 as described in the Introduction. The 
dotted curve belongs to S(a) obtained as the maximum 
taken over all occurring back-to-back circle pairs and all 
shifts p. Even for very large radii a, the background is not 
far below 0.3. Since there are 182 matched circle pairs, ev- 
ery peak in S(a) is generated by the circle pair with the 
highest signal. The figure also reveals the variation in the 
CITS signal for the individual pairs. For each pair the in- 
dividual CITS value is depicted by a small circle in figure 
[1] For q < 20° the CITS values of matched circle pairs 
are below the background of non-matched circle pairs, so 
that such small matched circle pairs can probably not be 
detected. Above a ~ 25° the matched circle pairs can be 
detected under the optimal conditions assumed here, where 
a pure CMB map without foregrounds and detector noise is 
available. 

The figure [T] also shows the CITS correlation S(a) ob- 
tained from the subset of circle pairs having all the same 
circle centres. For such a subset, the background due to 
non-matched circle pairs is much lower and not visible in 
figure [T] The CITS values of three subsets of such circle 
pairs are plotted with radii a = 20.5°, 24.7°, and 56.2°. The 
CITS values of the matched circle pairs (red dots in figure 
[TJ with a = 20.5° and 24.7° are below the full background 
(dotted curve). But in their neighbourhood towards smaller 
radii, there are correlations above the background niveau. 
This is a numerical observation which is confirmed in many 
instances: The maximal values in the CITS correlation are 
not obtained exactly on the matched circle pair, but, in- 
stead, for pairs shifted towards smaller radii. This implies 



that if one would discover a CITS signal in the real sky 
map, the radius a of the true circle pair is probably slightly 
larger than that corresponding to the peak. Along the axis 
of the selected 24.7° matched circle pair are three further 
pairs at radii 47.1°, 63.0°, and 76.9°, see the full curve and 
the corresponding red dots in figure [T] One observes that 
the shift towards smaller ring radii is less pronounced, if the 
ring radius a increases. 



3 SMOOTHING AND THE CITS AMPLITUDE 

Let us now turn to the important question how the smooth- 
ing of a CMB map can improve the amplitude of a CITS 
signal of a matched circle pair so that its peak emerges out 
of the background. This can be essential for detecting a pos- 
sible topology. There are several reasons why a smoothing 
of the map should be carried out. On the one hand, the tem- 
perature values of the pixels contain not only the pure CMB 
signal, but in addition the noise of the detector. Assuming 
that the noise is independent from pixel to pixel, a smooth- 
ing suppresses the noise contribution leading to a cleaner 
CMB signal. Even without noise a smoothing can be advan- 
tageous. Since the data are given on a discretised map, the 
scan strategy can miss the very localised peaks if the search 
grid is too coarse. If the data are sufficiently smoothed, even 
a relatively coarse search grid can find the topological sig- 
nal, and a lot of computer time can be saved. However, a 
too strong smoothing will eliminate the CITS signal. Thus 
there is a relation of how strong the smoothing should be 
for a given HEALPix resolution Ar s ;d e . As discussed in the 
Introduction, the CMB signal is composed from different 
contributions and only the Sachs- Wolfe contribution leads 
to a clear CITS signal while the Doppler and integrated 
Sachs- Wolfe contributions deteriorate it. The Sachs- Wolfe 
contribution provides the largest signal on scales around the 
first acoustic peak having a scale of roughly 0.8°. Thus the 
smoothing should at least be smaller than that in order to 
preserve the information due to the Sachs- Wolfe contribu- 
tion. 

We use a Gaussian smoothing which multiplies the 
expansion coefficients ai m of the temperature fluctuations 
5T(n) with respect to spherical harmonics Yi m (h) by 



(Hi. 



with 



7T 6>/180° 



(3) 



where 9 is the smoothing in degrees. 

This Gaussian smoothing is applied to the 3-torus sim- 
ulation with L = 1.5. The simulated CMB map is smoothed 
using the map in the HEALPix resolution A^idc = 4096. 
After the Gaussian smoothing, the map is downgraded to 
A^idc = 256, and the CITS correlation S(a) is calculated 
using this A^idc = 256 map. The result is shown in figure 
where the smoothing parameters 9 — 0.1°, 9 = 0.4°, 
and 9 = 0.6° are used. The panel (a) displays the corre- 
lation obtained from the 9 = 0.1° and 9 — 0.4° maps. The 
background of S(a) increases by increasing the smoothing 
parameter 9, but on the other hand, also the CITS sig- 
nal of the matched circle pair increases. In order to de- 
cide, whether there is an improvement by increasing 9 from 
9 — 0.1° to 9 = 0.4°, the figure also shows the difference 
5*6=0.4° (a) — S$=o.i°(oe) as a dashed curve. The difference 
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Figure 2. The CITS correlation S(a) is calculated for a 3-torus 
simulation with L = 1.5 using the HEALPix resolution N s ^c = 
256 and three different Gaussian smoothing. Panel (a) compares 
the CITS correlation for the two smoothing parameters = 0.1° 
and 6 = 0.4°, while panel (b) compares 9 = 0.4° and 6 = 0.6°. 
The dashed curve, which starts for a = at zero, displays the 
difference between the two CITS correlations S(a) shown in the 
panel. 



curve reveals the increasing background but also that the 
peaks are more pronounced for 9 = 0.4° , since the peaks are 
higher than the background. Thus the height of the CITS 
peaks increases more than the background. Increasing the 
smoothing further to 6 = 0.6° does not lead to an improve- 
ment as shown in panel (b). In panel (b) the difference curve 
Se=o.6° (a) — Sg=o. 4° (a) shows for ring radii a < 50° dips to- 
wards negative values. Thus the increase of the background 
is larger than the gain in the peak heights. Even for very 
large rings there is no genuine gain in the height of the 
peaks, but, instead, they are only extended towards smaller 
values of a. The focus of the smoothing should be put on 
the improvement of the signal for medium rings. Therefore, 
the smoothing parameter 9 should be smaller than 0.6°. We 




10 20 30 40 50 60 70 80 90 a 




10 20 30 40 50 60 70 80 90 a 

Figure 3. The CITS correlation S(a) is shown for the same 
simulation as in figure [2] but now calculated from the map in 
the HEALPix resolution N a ^ e = 512. Panel (a) compares the 
CITS correlation for the two smoothing parameters 9 = 0.1° and 
9 = 0.4°, while panel (b) compares 6 = 0.4° and 9 = 0.6°. The 
dashed curve, which starts for a = at zero, displays the differ- 
ence between the two CITS correlations S(a) shown in the panel. 



have investigated other smoothing parameters 9 in this way 
and conclude that for AWe = 256, the best smoothing pa- 
rameter is about 9 = 0.4°. 

Let us now address the best choice of the smoothing 
parameter 9 for the HEALPix resolution A^idc = 512. To 
that aim, the simulated CMB map is again smoothed using 
the map in the HEALPix resolution N s id e = 4096. But af- 
ter the Gaussian smoothing, the map is only downgraded to 
A^sidc = 512. The CITS correlation S(a) is calculated from 
this A^idc = 512 map and shown in figure[3]in the same way 
as figure [2] displays the A^ide = 256 results. Panel (a) shows 
the increase in the background by increasing the smooth- 
ing from 9 = 0.1° to = 0.4°, but at least for a > 40° 
the difference curve reveals an improved signal. However, 
for smaller ring radii the result is inconclusive since positive 
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as well as negative deviations from the background are ob- 
served. This suggests a twofold search if one insists on the 
-/Vsidc = 512 resolution: choosing 9 = 0.4° for larger rings 
and a correspondingly smaller value of 9 for smaller rings. 
The panel (b) shows that a further increase in 9 towards 
9 = 0.6° does not improve the signal as it was already found 
in the A 8x d e — 256 case. 

As a concluding remark of this section, it should be 
stated that a resolution of iV 8 ide = 128 should only be used 
to find matched circle pairs with large radii above a > 40°. 
A CITS search should thus be carried out in the A s id e = 256 
resolution using a Gaussian smoothing parameter 9 = 0.4°. 
If the sky map possesses no appreciable noise such that a 
downgrade to A s id c = 256 is not necessary as it might be 
the case for the Planck CMB map, the A s id c = 512 resolution 
would be preferable. 



4 THE SEARCH GRID 

In the CITS search every pixel of the HEALPix map is 
considered as a possible circle centre of a matched circle 
pair. The A s id e resolution leads to 12A s 2 idc pixels and de- 
termines how thoroughly the circles are scanned. However, 
it is not necessary to use for the circle centres, for which 
the CITS amplitude is to be computed, the same resolu- 
tion as for the CMB map. IVaudrevange et"al] <|2012T ) carry 
out their CITS search on a CMB map having the resolution 
A s idc = 512, but for the circle centres they use the coor- 
dinates of the pixel centres of the HEALPix map with the 
resolution A gr id = 128. This procedure saves a lot of com- 
puter time since the number of circle centres is reduced by a 
factor of 16. Except in this section, we always use the same 
resolution for the search grid and the CMB map. 

The coarser search grid harbours the risk that some 
matched circle pairs get missed. In order to test this issue, 
we calculate the CITS correlation S(a) from our 3-torus 
simulation with the HEALPix resolution A s id e = 512 and 
a Gaussian smoothing 9 — 0.4°. In the first case we use 
Ag r id = 512 and search all pixel centres of the A s id c = 512 
map, and in the second case only those of the A gr id = 128 
mesh. The result can be inferred from figure 2] where the 
dotted curve shows the A gr id = 512 search and the full 
curve uses the sparse A gr id = 128 grid. Several smaller rings 
are found only on the A gr id = 512 grid, see a = 24.7° 
in figure [4] The CITS correlations for the rings with radii 
a — 31.8° and a = 50.0° are reduced by a factor of two. 
Taking into account that our map contains the clean CMB 
signal without noise and residual foregrounds, it seems that 
also these signatures can be overlooked. Furthermore, if the 
multi-connected space has significantly fewer matched cir- 
cle pairs as it is the case for the L — 4 torus topology (six 
at a ~ 31° and three at a ~ 53°) mentioned in the In- 
troduction, even this back-to-back signature can be missed. 
To emphasise the large number of matched circle pairs in 
the L = 1.5 case, the figure [4] also displays the CITS val- 
ues at the pixel centres that provide the best match to the 
circle centres as small disks (A gr jd = 512) and open squares 
(A gr id = 128). It is obvious that a lot of matched circle pairs 
for the A gr id = 128 search are below the background even 
for a < 50°. All the open squares below the background 
level would never be considered as candidates for a possi- 
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Figure 4. The CITS correlation S(a) is calculated for the 3- 
torus simulation with L = 1.5 using the HEALPix resolution 
^side = 512 and a Gaussian smoothing = 0.4°. The difference 
in the two curves arises from the resolution Af gr id of the search 
grid which provides the coordinates of the circle centres. The 



dashed curve uses N, 



grid = 512 and the full curve N t 



grid 



128. 



Note, that both curves are computed from the same -/V si d c = 512 
CMB map. The values for the individual circle pairs are plotted 
as small disks (Af gr id = 512) and open squares (JVgrid = 128). 



ble refined search. Furth ermore, since the CITS search of 
IVaudrevange et "all (|2012t l refers to non back-to-back circle 
pairs, their background level is higher than in our back-to- 
back search. This worsen the situation, so that we doubt 
that their analysis exclude such matched circle pairs. 

Although the computer time would rise dramatically, 
one could ask how a search grid, which is finer than the 
resolution of the CMB map, would improve the CITS sig- 
nature. To address this point, the CITS correlation S(a) 
is computed for two search grids with A gr id = 512 and 
Agrid = 1024 on the 3-torus simulation with the HEALPix 
resolution A a id c = 512 and a Gaussian smoothing 9 = 0.2°. 
As figure [5] reveals the improvement justifies the increase in 
computer time only for small rings with radii a — 20° . . . 30° , 
where the stronger CITS signal is crucial. A similar compar- 
ison is carried out for a A s id c = 256 CMB map, where two 
search grids with A gr id = 256 and A gr id = 512 are used. It 
turns out that only for smoothing parameters 9 which are 
so small that the pixelization determines the resolution of 
the CMB map, an oversampling with A gr id > A s id c leads 
to significantly better results especially for ring radii in the 
range a = 20° . . .30°. 



5 NOISE AND THE CITS AMPLITUDE 

In the previous sections the CITS signature is analysed only 
for CMB map simulations which contain the pure CMB sig- 
nal. Only deteriorating effects that have their origin in the 
CMB physics such as the Doppler and the ISW effect were 
taken into account. The experiments, in contrast, yield CMB 
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Figure 5. The CITS correlation S(a) is calculated for the 3- 
torus simulation with L = 1.5 using the HEALPix resolution 
Aside = 512 and a Gaussian smoothing 8 = 0.2°. The figure shows 
how the CITS signal improves if a search grid with -/Vg r j<j = 1024 
is used on a N si(io = 512 CMB map. The results for both search 
grids are shown (full curve for Af gr id = 512 and dotted curve 
for iVg r id = 1024). The values for the individual circle pairs are 
plotted as small disks (JVgjjd = 512) and open squares (-/V gr id = 
1024). The difference curve is also shown, which fluctuates around 
zero. 



maps which are perturbed by noise and the beam profile of 
the detector. 

In order to address that issue, the 3-torus CMB map 
is modified according to t he beam prof i le by i,wi(l) of 
the Wl and VI channels ((Bennett et al . 2012) as given 
on the LAMBDA website using the HEALPix resolution 
iVgide = 4096. Thereafter, a downgrade from Aside = 4096 to 
A'aide = 512 is carried out. In the next step the noise of the 
pixel i is computed from the number A/ b 8 (i) of observations 
of the pixel i as a Gaussian random error with standard devi- 
ation o"o / v/ iVobs(*) and added to the map in the Aside = 512 
resolution. This map is expanded with respect to a; m and 
the Gaussian smoothing is realised as 

- l ( l + 1 )" 2 

aim = ai m e 2 / b V i,wi(l) , (4) 

so that the final smoothing is independent of the non- 
Gaussian beam profile b(l). This sequence of operations leads 
to CMB maps for the 3-torus topology comparable to the 
quality of the corresponding WMAP 9yr observations in the 
V and W bands. 

In this section CMB maps for the 3-torus with side 
length L — 4.0 are analysed. In contrast to the L = 1.5 
simulation, only multipoles I with I ^ 1000 are included in 
the L = 4.0 case. This simulation is accurate enough for an 
analysis with respect to noise as can be inferred from figure 
[6] For I > 1000 the angular power spectrum 5Tf is noise 
dominated for the V and W bands of the WMAP measure- 
ments which are put into focus in this section. The compari- 
son of the two panels of figure [6] shows that the W band has 
a higher noise contribution than the V band, although the 




500 1000 I 1500 



Figure 6. The angular power spectrum <5T ; 2 = 1(1 + l)C;/(27r) is 
shown for two 3-torus simulations with L = 1.5 and L = 4.0. Two 
power spectra belonging to L = 1.5 and L = 4.0 take only the 
beam profile into account. These two curves are nearly indistin- 
guishable and drop towards zero for large values of I. A marginal 
difference is visible in the W band simulation for I > 1000, where 
only the L = 1.5 simulation has power. Furthermore, the curve 
increasing for large / belongs to the L = 4.0 simulation where in 
addition to the beam profile also the detector noise is accounted 
for. 

W band beam profile reveals tentatively the third acoustic 
peak by ignoring the noise. 

It is interesting to evaluate the CMB maps, which are 
processed to mimic the V and W band maps of WMAP, 
in the AWe = 512 resolution without further smoothing. 
Figure [7] shows that the CITS signature of the 6 ring pairs 
at a = 31.1° is missing at the A^ide = 512 resolution in 
both bands. Only the 3 ring pairs at a = 52.7° emerge 
slightly from the background, where the signature is more 
pronounced in the V band having less noise. The CITS sig- 
nature is swamped by the noise as a simple downgrade to 
the A s idc = 256 resolution reveals. Here the average over 4 
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Figure 7. The CITS correlation S(a) is calculated for the 3-torus 
simulation with L = 4.0 using the beam profile and the detector 
noise of the V band in panel (a) and W band in panel (b). The 
iVgide = 256 map is obtained by a downgrade of the -/V s j<j c = 512 
map. This topology has 6 ring pairs at a = 31.1° and 3 ring pairs 
at a = 52.7°. 



Figure 8. The CITS correlation S(a) is calculated for the 3- 
torus simulation with L = 4.0 using the beam profile and the 
detector noise of the V band. The CMB map is smoothed using 
the parameters 6 = 0.4° and 6 = 0.5°. In panel (a) S(a) is 
obtained from the N s ^ c = 512 CMB map, while panel (b) shows 
the corresponding curves obtained from the N^c = 256 CMB 
map. The difference curve is also shown. 



pixels yields the temperature value in the N s id c = 256 reso- 
lution. The improvement can be seen by the dotted curves 
in figure [7] Despite the increased background, the CITS sig- 
nature is visible more clearly. This demonstrates that the V 
and W band CMB maps have to be smoothed in order to 
reduce the degrading effects due to noise. 

The importance of the smoothing to reduce the noise 
is also revealed by the following analysis. At first the CMB 
map with the V band properties in the iV s id G = 512 reso- 
lution, whose CITS correlation is shown in figure [7{ a), is 
smoothed according to eq. ((4]). This smoothing is done in 
the iVaido = 512 resolution for various values of 9. The opti- 
mum is found in the range between 9 — 0.4° and 9 = 0.5°. 
For smoothing parameters 9 below 9 = 0.4° the noise de- 
grades the CITS signal and above 9 = 0.5°, the smooth- 



ing removes too much CMB information. The CITS corre- 
lation is plotted in figure [8j a) for 9 — 0.4° and 9 = 0.5°, 
and the improvement compared to figure [TJa) (full curve) 
is obvious. In addition, the difference curve is plotted which 
shows that both smoothing parameters lead to nearly indis- 
tinguishable results. In panel (b) of figure [8] the analogous 
curves are shown for the CMB map which is downgraded 
to the JVsidc = 256 resolution. Since the computer time is 
reduced for the smaller iV s ide — 256 resolution, we will use 
in the following this value of N s id e . 

Up to now the analysis is based on a single CMB sim- 
ulation for the 3-torus topology either for L — 1.5 or for 
L = 4.0. However, the CMB simulations require Gaussian 
random initial conditions leading to a statistical ensemble 
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Figure 9. The CITS correlation S(a) is calculated for 100 differ- 
ent 3-torus simulations with L = 4.0 using the beam profile and 
the detector noise of the V band. The CMB maps are smoothed 
using the parameters 9 = 0.4° (panel (a)) and 8 = 0.5° (panel 
(b)) and downgraded to N sidc = 256. The blue dotted curve shows 
the maximum of S(a) taken over all 100 simulations. The red full 
curve shows the corresponding minimum and represents thus the 
worst case occurring in the 100 simulations. The shaded band en- 
compasses all values of 5b g (a) occurring in the 100 simulations, 
while the grey band contains 66 percent of the Sb g (o) values 
around their median. This reveals the variation of S(a) with re- 
spect to the random initial conditions of the CMB simulations. 

of CMB simulations belonging to a given topology of a fixed 
side length L. In order to address the question how this vari- 
ability affects the CITS signature, an ensemble of 100 CMB 
simulations for the 3-torus with L = 4.0 is generated. This 
leads to 100 CMB maps where the beam profile bvi{l) and 
the detector noise of the V band is taken into account as de- 
scribed above. The results for the two smoothing parameters 
6 = 0.4° and 6 = 0.5° are presented in figure H The dotted 
curve shows the largest CITS signature that occurs in the 
100 simulations, while the full curve shows the worst case 
that is the minimum over the 100 CITS correlations S(a). 
Individual matched circle pairs can possess even lower val- 
ues, but note that S(a) takes the maximum of the 6 pairs at 
a = 31.1° and of the 3 pairs at a = 52.7°. A detailed anal- 
ysis reveals that the smoothing parameter 9 = 0.5° yields 
marginally better peaks than the choice 6 = 0.4°. 

At radii a where no matched circle pair is nearby, the 



correlations Sij(a,f3), see eq. (JTJ, show (positive) correla- 
tions and (negative) anti-correlations with the same ampli- 
tude, so that the minimum of the anti-correlations, defined 
by 

S be (a) = max ( - ^(o.jS)) , (5) 

is a measure of the background even at radii a where 
matched circle pairs do occur. Thus it quantifies the detec- 
tion threshold. The figure [5] shows the distribution of Sb s (ct) 
obtained from the 100 simulations as the shaded band. This 
distribution is asymmetric with a tail towards larger val- 
ues, and the figure [5] also shows as the grey band the band 
which contains 66 percent of the values of S'bg(a) around 
their median. 

One observes in figure [U] that the CITS signature of 
the six matched circle pairs at a — 31.1° is only slightly 
above this background in the worst case simulation. Since 
the measurements contain residual foregrounds not taken 
into account by our simulations, which further deteriorate 
the signal, it seems that even the detection of matched circle 
pairs around a = 30° cannot be assured. For larger ring radii 
a, see the signature at a = 52.7° in figure [9l the situation is 
far better, and such matched circle pairs should be detected. 



6 CITS SIGNAL FOR THE W AND V BANDS 
USING MASKS 

The noise is not the only problem for the CITS signature, 
since foreground sources dominate the CMB in some parts 
of the sky so that the CMB signal cannot be reliably re- 
constructed. The way out found in the literature consists of 
masking the problematic sky regions in the CMB map and 
replacing the temperature values by those of the ILC map. 
Here, we do not use the pixels which are too severely con- 
taminated by foregrounds and compute the CITS correlation 
only with those pixels which are lying outside the KQ75 and 
KQ8 5 masks provided by the WMAP team (|Bennett et al.l 
|2012J). Thereby arises the problem that the values of the 
CITS correlation cannot be compared for circles having the 
same radius a, since the correlations are computed from a 
different number of pixels depending on the orientation of 
the circles with respect to the mask. In a given HEALPix 
resolution JV a id e , a circ le of radius a has t he ring index n r i ng 
in the RING scheme (|G6rski et al.ll2005T l of the HEALPix 
map, i.e. firing ~ n T i ns (a). Define the number n e ff as the 
number of pixels that are used for the computation of the 
CITS correlation, that is n e g n^f* = 4min(n r i n g, N s id e ) 
for n r i n g ^ 2iV s idc with equality for the case that the mask 
does not affect any pixel on the ring. The correlation S(a) is 
then analysed on the two-dimensional grid n e g versus n r i ng . 

According to the results of the precedin g sections, the 
foreg round reduced V and W band maps IjBennett et al.l 
120121 ) are expanded in the resolution iV s idc = 512 with re- 
spect to ai m and then are smoothed according to eq. (Q. 
This yields the CMB maps which are analysed in the fol- 
lowing subjected to the KQ75 or KQ85 mask. 

We analyse the CMB maps with the resolution N s id c = 
256 for the smoothing parameters 6 — 0.4° and 6 = 0.5° by 
computing the CITS correlation S as a function of n r i ng and 
n e s- We do not find a pronounced peak in the data which 
could be interpreted as a signal for a matched circle pair. 
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Figure 10. The CITS correlation S is computed for the V and 
W bands of the WMAP data using the resolution -/V s idc = 128 as 
described in the text. The maps are smoothed to 0.5° and sub- 
jected to the KQ75 mask. The value of S encoded by brightness 
is shown as a function of n r ; ng and n e fj. Two interesting values 
are marked by arrows. 



Although our tests with the topological simulations dis- 
cussed in the previous sections suggests otherwise, we nev- 
ertheless searched in the lower resolution AT a jde = 128 after 
peaks. This differs from the search in the Af 8 ide = 256 maps 
by putting the average of four AT B ide = 256 pixels into one 
A^sidc = 128 pixel. By this additional averaging the detector 
noise is suppressed more than necessary according to our 
simulations. However, a possible residual foreground signal 
on small scales would also be smoothed out thereby reduc- 
ing its contribution to the m-weighted circle signature de- 
fined in eq. HJ. The figures \W\ and [11] present the results 
for the CITS search in the N s idc = 128 maps. In figure \W\ 
the value of the CITS correlation is encoded by brightness. 
As already discussed, the horizontal axis n r in g corresponds 
to the radius a of the circle in the iV S ide = 128 resolution, 
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Figure 11. The CITS correlation S(a) shown in figure 1101 is 
shown as a one-dimensional plot by taking the maximum of S 
over the stripe n e g £ [175, 250] in panel (a) and over the stripe 
n c ff £ [250,500] in panel (b). The peak at n r ; ng = 202 visible in 
panel (a) is marked by the right arrow in figure lTOl The left arrow 
corresponds to the peak at n r ; ng = 72 seen in panel (b). 



while the vertical axis over n e g separates the CITS values 
with respect to the number of pixels on which its calcula- 
tion is based. Towards smaller values of n e g, the CITS value 
and correspondingly the brightness increases, in general, be- 
cause a lower number of pixels enhances the probability of 
a large correlation by chance. A pronounced peak occurs at 
Wring = 202 corresponding to a — 73.7° and n e g = 188. 
Its position is marked by an arrow in figure [10] A further 
peak at the position n r i ng = 72 corresponding to a smaller 
ring with a = 26.6° and n er r = 256 is also marked. Al- 
though CITS values belonging to different values of n e ff can- 
not be strictly compared due to their different significance, 
we nevertheless combine two stripes with n e fi £ [175, 250] 
and n e g € [250, 500] by taking the maximum in the given 
n e ff interval for a fixed value of firing- The results are plot- 
ted in figure[TT]for both stripes for the V and W band maps 
subjected to the KQ75 and KQ85 masks. Panel (a) reveals 
the pronounced peak at n r i ng = 202. The signal is visible in 
all 4 cases shown in the figure. A further peak in panel (a) 
at n r in g ~ 90 emerges only by using the KQ75 mask and is 
absent using the KQ85 mask. The decrease in the CITS cor- 
relation is thus caused by pixels close to the boundary of the 
KQ75 mask. This case demonstrates that a few additional 
pixels can destroy a CITS signal. The peak at n r i ng = 72 
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Figure 13. The full curve shows the CITS correlation S obtained 
from the addition of the V and W band CMB maps, while the 
dotted curve belongs to a torus simulation having a side length 
L = 4.0. The maximum of S over n e g ^ 300 is plotted. The reso- 
lution N a idc = 256 together with the KQ85 mask and a smoothing 
with 9 = 0.5° is used. 



Figure 12. The W band temperature fluctuations subjected to 
the KQ75 mask are shown together with the position of the circle 
pair with n r ing = 72 (radius a = 26.6°) in panel (a) and with 
Wring = 202 (radius a = 73.7°) in panel (b). 

shown in panel (b) is more clearly seen in the W band data 
and is independent of the mask. 

It is instructive to discuss the position of the circle pair 
on the CMB map. The large circle pair with n r i ng = 202 lies 
in the neighbourhood of the galactic plane as figure 112T b) 
shows. This is the reason why large parts of this circle pair 
are excluded by the mask. In addition, this location bears 
the risk that many pixels used in the computation of the 
CITS correlation contain significant amounts of foreground 
contributions. The situation is more favourable in the case of 
the smaller circle pair with n r i ng = 72 shown in figure [127 a). 
This pair is far away from the galactic plane and most pixels 
are taken into account. At n r j ng = 72 one has n^f x = 288 
which is close to n c ft = 256 found for this circle pair. On 
the other hand, the simulations of the preceding sections 
demonstrate that a circle pair with such a small radius is at 
the detection threshold. 

The CITS signatures in the resolution -/V s id c = 128 
at n r i ng = 72 and n r i ng = 202 are probably generated by 
chance. This is suggested by a comparison with the results 
in the resolution A s id c = 256 where these circle pairs reveal 
no unusual high correlation values. Furthermore, we add the 
V and W band maps after correcting for their different beam 
profiles. After smoothing the summation map with 9 = 0.5°, 
it is downgraded to A s id e = 256. The maximum of the cor- 
relation S with n e g 300 is plotted in figure [13] as the full 
curve, where the KQ85 mask is used. The values of n r i ng for 
iV s idc = 256 are twice that of A S ide — 128 for the same radii 
a. So peaks around Tiring = 144 and n r i ng = 404 would be 
expected from the A s idc = 128 results. Although a peak at 
Wring = 140 appears, it has nothing to do with that in the 
Aside = 128 resolution, since its circle pair has a different 
orientation. Furthermore, this peak is absent if the KQ75 



mask is used. The figure [13] also shows as a dotted curve 
the CITS correlation for a torus model with L = 4.0 having 
circles at a = 31.1° (n r m g = 168 for A s id c = 256) and at 
a = 52.7° (nring = 279 for A sidc = 256). The large circles 
at n r i ng = 279 are clearly revealed, while none of the six 
small circle pairs at n r i ng = 168 leads to enhanced correla- 
tion values. The fact that the circle pairs at a = 31.1° are 
not detected, suggests that the peak at n r i ng = 140 corre- 
sponding to a = 25.8° is produced by chance. Furthermore, 
this fact also shows that matched circle pairs should have 
radii larger than a — 30° in order to be detectable with 
WMAP data. 



7 SUMMARY AND DISCUSSION 

The circles-in-the-sky (CITS) signature provides a test for 
a possible non-trivial topology, i. e. for a multi-connected 
space, of our Universe. The detection threshold of a matched 
circle pair depends on the radius a of the circle. In order to 
determine the detection threshold, accurate CMB maps are 
simulated for the 3-torus topology, and their CITS correla- 
tions are computed for various circumstances. 

The CITS correlation is introduced in eqs. ^ and ([2]| 
in section [5] This topological detection measure is applied 
to a simulated 3-torus CMB map for a restricted search for 
back-to-back circle pairs. It turns out that for radii a > 25° 
the CITS signal is strong enough for a detection of matched 
circle pairs. Furthermore, it is found that the highest cor- 
relation values are not found on the matched circles, but 
rather for neighbouring circles having slightly smaller radii. 
This effect is the more pronounced, the smaller the radii of 
the matched circle pairs are. 

The experiments provide CMB maps in which the res- 
olution is determined by the size of the pixels. Section [3] in- 
vestigates the question how an additional Gaussian smooth- 
ing of the CMB data can enhance the detection probabil- 
ity for a matched circle pair. The amplitude of the CITS 
signal is compared with the detection threshold for several 
CMB simulations with different HEALPix resolutions Aside 
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and different Gaussian smoothing parameters 9, see eq. @. 
This analysis leads to the conclusion that for A^ide = 256 
a smoothing parameter of 9 = 0.4° is preferable if the pix- 
els contain the clean CMB signal. But for A^ide = 512 one 
should prefer a smaller smoothing parameter 9 in order to 
detect also matched circle pairs with smaller radii. Experi- 
mentally, such clean sky maps are not available because of 
detector noise and residual foregrounds. 

Before these problems are addressed, we turn in section 
[4] to the question how sharp the peaks are which belong to 
genuine matched circle pairs. This determines the resolu- 
tion of the search grid AT gr id whose coordinates are used as 
possible circle centres for which the CITS correlation is com- 
puted. It is demonstrated that for a smoothing of 9 = 0.4° 
several matched circle pairs get missed at iVg r id = 128 while 
using iVsidc = 512 for the CMB map. Although one can 
find matched circle pairs at iVgrid = 128, one cannot exclude 
matched circle pairs if no circle pairs are found at this search 
grid resolution. 

In sections [2] to [4] only general properties are studied 
without reference to a specific experiment. In section [5l 
however, CMB maps for the 3-torus topology are generated 
which take the noise properties of the V and W band chan- 
nels of the WMAP satellite into account. The CITS correla- 
tion for a CMB map having the quality corresponding to the 
WMAP 9yr observations is computed in the Af s id c = 512 res- 
olution. It turns out that CITS signal of the matched circle 
pairs is either absent or very faint which demonstrates that 
the noise is too severe in order to search directly in the orig- 
inal V and W band maps for a topological signature. This 
result is further substantiated by smoothing the A^ide = 512 
map using the parameters 9 = 0.4° and 9 — 0.5°. These 
additionally smoothed AT S ide = 512 maps reveal the CITS 
peaks at the positions of the matched circle pairs. Such a 
smoothing is necessary for a sufficient suppression of noise. 
Since the above conclusions are based on a single CMB sim- 
ulation, a set of 100 CMB maps is generated for the 3-torus 
topology with side length L = 4. Their 100 CITS corre- 
lations reveal the variability with respect to the Gaussian 
random initial conditions which determine how the eigen- 
modes are weighted. As shown in figure |9j the variability of 
the height of the CITS peaks at circle radii around a ~ 30° 
is so important that some matched circle pairs are at the de- 
tection limit. Around a ~ 50° the variability is of the same 
order, but the detection probability is not affected due to 
the more pronounced peaks. 

Section [6] addresses as a further complication the resid- 
ual foregrounds which can be so severe that it is advisable 
to mask the unsafe regions. As a consequence the CITS 
correlation can only be computed from the unmasked frac- 
tion of the circle leading to a statistical significance de- 
pending on the number n e ft of pixels used for the compu- 
tation of the CITS correlation. The search for a CITS sig- 
nature in the foreground reduced V and W band maps with 
N sidc = 256 and 9 = 0.4° or 9 = 0.5° subjected to the 
KQ75 or KQ85 mask yields a negative result. Although the 
results obtained from our simulation would not advise to 
use the A^ide = 128 resolution in a CITS search, we nev- 
ertheless carry out the Aside = 128 search and find some 
interesting peaks. A peak emerges from the background at 
the large circle radius a = 73.7° as seen in figures [10] and 
111( a). The corresponding circle pair lies, however, as seen 



in figure [T27 b). very close to the galactic plane. Since one 
expects the largest residual foregrounds close to the galactic 
plane, the suspicion arises that the signal might be spu- 
riously generated by residual foregrounds. However, if one 
accepts such an argument, the reverse argument that these 
very residual foregrounds can also destroy a true CITS sig- 
nal, have to be accepted in the same way. Furthermore, since 
so many pixels are masked in this circle pair, the high CITS 
correlation can also be produced by chance. A second inter- 
esting peak in the CITS correlation belongs to a circle ra- 
dius a — 26.6°, which is also probably generated by chance 
as our simulations show. On the other hand, this circle pair 
is far from the galactic plane, and few pixels are masked 
and the residual foregrounds are less important. As demon- 
strated in section [6l the six circle pairs at a = 31.1° in the 
L = 4.0 torus topology are not detected by using the quality 
of the WMAP data, which shows that matched circle pairs 
have to be larger than at least a = 30° in order to be de- 
tectable. We are thus eagerly awaiting new more accurate 
data in order to see, whether the CITS peaks will with- 
stand or prove to be a spurious result. Even if these CITS 
signatures will not endure, a multi-connected space can pos- 
sess circle pairs that are not back-to-back and these are not 
considered in this paper. Examples for topologies with such 
generic m atched circle pairs are p rovided by the half-turn 
space E2 (jAurich and Lu stie 2011) a nd by s pherical double- 
action manifolds ( Aurich and Lustidl2012b1 lah. These spaces 
would require an extended search for generic matched circle 
pairs. 
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